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ijr^ ' In this paper we develop a formalism to describe a superfluid in a gravitational 



background. This formalism is based on a covariant generalization of the field 
description for a superconductor in terms of a U(l) spontaneous symmetry 
breaking. We study the stability of the solutions for a vortexless fluid and 
the force acting on vortices in the fluid, which is a generalization of the well- 
known flat space-time Magnus force. To clarify the development we include 



^ ' the explicit discussion of two particular cases, one of them of astrophysical 



interest. 
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I. INTRODUCTION 

In this paper we study the behavior of a superfluid in presence of a gravitational field. 
Under the usual conditions that can be obtained in a laboratory only helium shows the 
superfluid effect, because this is the only substance that remains liquid at temperatures low 
enough for quantum effects to dominate. New and appealing scenarios with very interesting 
possibilities are provided by some astrophysical systems. One scenario is realted to neutron 
stars, where there are very high pressures and densities, and thus where superfluid effects can 
appear at relatively high temperatures. A widely accepted description for the constituent 
matter of a neutron star considers two coexisting quantum superfiuids, neutron and proton 
seas, where quantum dynamics plays an essential role |l|. But, in contrast with laboratories 
on Earth, there we have not only superfiuids but also strong gravitational fields. This makes 
it possible to have significant gravimagnetic effects, which could produce phenomena similar 
to the ones that appear in the case of the superconductors [0]. We can also mention a 
more speculative scenario, of interest at a cosmological level, which is the formation of a 
condensate in the sea of relic neutrinos 0. Scenarios like these open the possibility of 
gaining access to a phenomenology which would allow a deeper understanding of different 
aspects of the gravitational interaction in quantum systems. 

There have already been some attempts to describe a covariant superfluid, usually stress- 
ing classical geometrical aspects of the problem @] . We present here a different approach, 
which is a covariant generalization of the description of the superconductor phenomenology 
based on a U(l) spontaneous symmetry breaking p. This approach allows us to discuss 
dynamical aspects of the superfiuids in a gravitational background, on a very general and 
well established basis, with a clear phenomenological interpretation. 

In the following section we develop the general formalism for a superfluid in a gravita- 
tional background, and discuss the dynamics of the scalar bosonic excitations. We make 
explicit the corresponding equation of motion and the conditions that warrant the stability 
of the ground state. In Section 3 we consider the presence of vortices and study the forces 



acting on them, using a Kalb-Raniond description for the vortices [0. In particular, we 
obtain the force due to the background, which is a generahzation of the well known Magnus 
force, and the stable vortex configurations. Section 4 contains the analysis of two cases, 
a superfiuid in a laboratory on Earth and in a neutron star, which exemplify the general 
formalism. Finally, in the last section we make some comments and remarks. 

II. SUPERFLUID IN A GRAVITATIONAL BACKGROUND 

The superfluidity phenomena is related at the microscopic level to the formation of a 
condensate that spontaneously breaks a global symmetry. The spontaneous breaking of the 
symmetry leads to a Nambu-Goldstone excitation with zero energy in the limit of vanishing 
momentum. We develop here a presentation for the superfiuid that follows the one for 
superconductors in flat space-time given in Ref. [^]. We start with a fluid which is formed 
by particles represented by a matter field \1' with the simplest U(l) global symmetry group 
"^{x) — > e*^\l/(x). To be specific we can think of the neutron fluid in a neutron star and 
the baryonic symmetry that leads to the baryon number conservation in the system. We 
can always perform a field splitting 

^(x) = e*<^(^) tPix) . (1) 

The group transformation acting on the boson is (f){x) — > 0(x) + A. 

The U(l) invariant density Lagrangian must be a function of the derivatives of (^ and 
the U(l) invariant field ip. The most general Lagrangian density allowed by the symmetries 
is a nonlocal function of the field, where the nonlocality extends over a range of the order of 
the coherence length of the superfiuid. Given that we are interested in the macroscopic fluid 
motion, we will only consider local terms in the Lagrangian that effectively describe the long 
range behavior. Such terms must be scalars and should be constructed as a contraction of 
covariant quantities. The only possible fields involved are the gradient of (p, the metric and a 
number of tensors A'^, A^'^, etc., which depend on the field ip and must satisfy the additional 



requirements of symmetry eventually imposed by this field. Thus, the Lagrangian can be 
written as the expansion 

£ = £4^] + \^[^lj]ct^,, + ]^\^'[4j] <\>^.<\>^, + ... , (2) 

where £<, is a independent Lagrangian. The symmetry breaking implies that the Goldstone 
field has an effective dynamics below some typical energy where the condensate forms and 
where the tensor \^^ develops a non null expectation value. The dynamics of the superfluid 
is completely described by the Goldstone boson because the gap makes i\} nondynamical at 
low energies. However, the expectation values of the tensors \^^ and A^ generally depend 
on external forces, such as the ones induced by the gravitational field and the ones on the 
background state. That is, in general \^^ and A'^ should satisfy classical equations that 
involve the gravitational field, of the type given by the hydrostatic equilibrium, while the 
Goldstone boson should satisfy the full quantum dynamics given by the Lagrangian (||). 
Besides this, if the Goldstone field satisfies a four dimensional dynamics the A^^ tensor 
must be non degenerate. The effective theory makes contact with the microscopic theory 
through the values of these tensor coefficients. Given a particular case these quantities 
could satisfy additional symmetries, which would allow ad hoc approximations to achieve a 
concrete physical characterization. 

The field is dimensionless and thus the tensor A'^i'^^-.-^n must have a (4 — n) dimension, 
being typically of the order of /i^"*^, where /i is the condensate scale. The low energy require- 
ment states that the typical scale e of (9j,0 satisfies e <^ /i, which justifies us in expressing 
the Lagrangian as an expansion on the derivatives of and in cutting this expansion at 
order two. In the case where the system contains several different energy scales a sensible 
expansion can involve different powers of the space and time derivatives [0. This could 
be the case for laboratory superfluids where the Fermi momentum and the fermion masses 
are very different. We will focus on superfluids such as the one in a neutron star, where 
these quantities are similar. For the description sketched above to be valid during the whole 
evolution of the system, the gradient of must not only be small at a given instant, but 



also at all times. From the equation of motion of 0, this implies that the coefficients A'^^ 
and A^ must be sufficiently smooth. 

The original theory is invariant under general coordinate transformations and redefi- 
nitions of the field factorization. Once the background has been established the effective 
theory admits only a restricted refactorization. The effective second order Lagrangian will 
be invariant under the field transformation 0^0' = — A,-?/'— ^ ■?/;' = e*^^^V > if ^,v is of 
the same order as (p^u- 

Inspired by the physical picture for the neutron star matter, where neutron pairs have 
nonvanishing expectation values, we will assume that the U(l) baryonic symmetry is broken 
to Z2, the subgroup of transformations with A = and A = vr. This is similar to the 
superconductivity effect, where there is a spontaneous breaking of the U(l) electromagnetic 
gauge symmetry to Z2. As parametrizes U(l)/Z2, and + vr are taken to be equivalent. 

The baryonic current is given by 

BT 

/ = ^ = A^^(x)0,, + A'^(x). (3) 

It contains two contributions. One is given by the Goldstone field, in terms of the derivatives 
of 0, and the other comes from the background baryonic current given by A^. The continuity 
equation for the baryonic current, d^{y/gj'^) = 0, is equivalent to the equation of motion of 
the Nambu-Goldstone field 

D^(A'^'^(x)0,, + A^(x))=O. (4) 

The velocity of charge transport can be defined by the relation j'^ = UoU^-, where Uo is 
interpreted as the charge density in the fiuid rest frame 0. The superfiuid state without 
vortices in a fiat space-time is characterized by the relation V x v = 0, where v is the 
velocity of the fiuid. In a curved space-time, this relation is replaced by the one obtained 
from taking the curl of the gradient in Eq. (|^) 

e^^^'Ds(\-^Y{noU^-\^) = Q. (5) 



This last equation, together with Eq. (H) defines the evolution of a vortexless superfluid in 
a gravitational background. 

The energy-momentum tensor is 



rpfl 



The density of force acting on the fiuid due to the interaction with the background is /^ = 
DyT^^ = —DI^C, where D^ means the derivative with respect to the explicit dependence of 
the Lagrangian in x'^. Thus 

We can define a global time by a time-like vector ^^ ~ dt, with ^.^ = 1, orthogonal to 
space-like surfaces that are set to be equal time three-dimensional subspaces. The density 
of energy in such equal time space-like subspaces is p^ = i^Tj^i": 

PC = ii. A^'^r - ^A'^- ) 0,^0,. + ((A -or- A'^) 0,. , (8) 

and the power exchanged with the background is given by the projection of f^ along ^^: 

w = ri, = ~D^{X'^) <P,^<P,r-D^i\'^) 0,., (9) 

which depends on the directional derivatives of X'^'^ and A'^ along ^. 

The covariant statement of the existence and stability of the ground state can be ex- 
pressed in terms of a minimum of the energy. If we want a local extremum for o- = 0, it 
must be 

(A-Oe. = A.. (10) 

This implies that ^o- must be proportional to Ao-, the baryonic current for 0.o- = 0. Thus the 
background explicitly breaks the covariance of the theory. Furthermore, if the extremum 
corresponds to a minimum, the matrix 

M"' = ^£^ = (ep Cg: + i, Cg: - gX) ^'^ W 
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must be positive definite. As stated above, if the Goldstone boson lias a true four dimensional 
dynamics, A^ is regular, and thus it has one time-like eigenvector and three space-like ones. 
If we also require that the background be invariant under the time reversal transformation 
C,cr -^ —^a , the M^^ matrix becomes reducible into two non null submatrices: a one- 
dimensional one in the subspace spanned by C,a^ 

i^M^X, = iA^X, , (12) 

and a three-dimensional one in the space-like subspace orthogonal to ^o- 

(e^e. - 9'.) M^^ (ci, - (?;) = - {ei. - €) A'^^' (re. - ^;) • (13) 

Therefore, in this case to have an energy minimum, the component C,a^'^^C,ri of X'^^ must be 
positive, whereas the matrix of the components in the orthogonal subspace must be negative 
definite. Thus the minimum energy condition guarantees dispersion relations with real prop- 
agation velocities for the low-energy excitations. If the space-like submatrix is completely 
degenerate, the propagation velocities are the same in all the spatial directions and we have 
an isotropic background. Otherwise we have an anisotropic one, with different propagation 
properties for the low energy excitations according to the different spatial eigenvectors of 

The spacial components Ui of the velocity are proportional to — Af 0^^ and this implies 
that the low energy expansion is also a small velocity expansion -Uj ^ 1. 

We can implement a hydrodynamical description of the fluid by writing the energy mo- 
mentum tensor in terms of the current. It reads 

Ti: = y {u^ {\'')^^ {n,u^ - A") + u, (A-^)^ (n^w" - A")) 

-gii^[X'^)^^{nyu^-X^X^) . (14) 

This tensor only contains the terms corresponding to the dynamics of the excitations with 
respect to the background. To have the energy momentum tensor of the fluid as a whole, it 
is necessary to add the contribution of the rest of the fluid. Assuming that the fluid is an 
isotropic perfect one up to first order in the velocities, we have 
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(P + P) ri, 



f (^-'); . (15) 



no 2 

where /x is (for zero temperature) the chemical potential of the global conserved number, p 
is the energy density and p is the pressure of the fluid. 

III. VORTICES AND THE GENERALIZED MAGNUS FORCE 

We have not considered the existence of vortices up to now, but they will be present if 
the fluid has some angular moment. In general, the superfiuid cannot rotate as a whole, 
and therefore if it has a non null angular momentum it must be supported by vortices. It is 
difficult to analyze the vortex-fiuid interaction in the Goldstone field representation for the 
superfiuid, because the vortex imposes non trivial boundary conditions to the scalar field. 
In the large distance limit where the vortex size could be considered negligible these con- 
ditions manifest themselves as a field nondifferentiability that leads to a noncommutativity 
of its partial derivatives. A more suitable approach could be developed relying on the dual 
description of the scalar field in terms of a Kalb-Ramond one, an antisymmetric tensor of 
second order B^^ ^^- Both massless theories, a k-form field and a (n-k-2)-form field, where 
n is the space-time dimension, are equivalent not only at a classical level, but also at the 



quantum one. This point will be discussed in detail elsewhere [|T0|- We will use a generalized 
version of this duality here. 

The dual variables are defined by expressing the Noether conserved baryonic current as 
a topologically conserved current. This is achieved with the identification 

jx ^ ]^E^-P^H^,p , (16) 

where H^y„ = d^B^^ + di^B^^ + d^jB^y and E is the Hodge tensor with covariant and 
contravariant components E^^px = \/—gSp,yp\ and E^^''''^ = —^^e^^^^ respectively. The 
Levi-Civita symbols denote the index permutation signature. The completely antisymmetric 
tensor E commutes with the covariant derivatives. The inverse relation to (0) is H^^yp = 

^fiupX'J ■ 



The divergence of the current J^ is proportional to the exterior derivative d A H, and 
nuUifies because H = d A B. In components it reads 

e'^'^p^dxH^^p = . (17) 

In this way the equation of motion for the field </>, which is equivalent to the baryonic 
current conservation, is satisfied automatically in a topological sense in the Kalb-Ramond 
field scheme. 

In a similar way, the topological property of the scalar field d A {d A(j)) = becomes the 
equations of motion for the Kalb-Ramond field -B^i/, by means of the duality transformations. 
Duality states a correspondence between the physical solutions of the two theories. This 
correspondence is realized as a canonical transformation at the level of Dirac brackets, which 
define the quantum theory |[10[| . 

Both descriptions are linked by Eq. (p^, which is a relation between the field derivatives 
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gE^^^'^i7,.,-A-JA;,i , (18) 

where the tensor A~^ is the inverse of A'^'^. This relation is a generalization of the one used in 
Ref . [0 , where only the particular case Ao-^ oc Qo-k is considered. The curl of ^ is null, and 
this implies that the curl of the right hand side must also be null. This gives the equations 
of motion for the B field: 

E-^^^D, i^^E^'^P'^H.^pX-l - KlX^'^ = . (19) 

The corresponding action that leads to this equation can be written as 

Sb = ^J d'x^/^ H^p,E^P'<^ Q^E^''P^H,,,\-l - KlX'^^ . (20) 

Equating the energy-momentum tensor from this action with the one for the scalar field, 
Eq. (^, using the relation (|18|), we get k = —\. 

The 2-form Kalb-Ramond field naturally couples with the two dimensional manifold 
generated by the world sheet of the vortex moving along time. This formalism is advanta- 
geous because it allows us to describe the fluid-vortex interaction with an explicit term in 
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the action, without using any extra boundary conditions. Using the antisymmetric tensor 
jtiu _ fj J dg-fj''^ §^ (^x — ylcr^r)), which describes the surface that sweeps the vortex in space 
time, with the coordinates (a, r) parametrizing the world sheet, the fluid-vortex interaction 
term is 

C = B,,J^'' , (21) 

a metric independent term because it corresponds to the surface integral of a two form. It 
is possible to think of different interaction terms that depend on the metric or that have 
field derivatives, but for dimensional reasons those terms should be suppressed. They must 
contain additional powers of the superfiuid coherence length over the vortex curvature scale 
or the space curvature, and can be neglected in this context. 

Taking into account the interaction with the vortex, the equation of motion is written 

^E^^^'^D, (^^E^'^'^^H.^.X-'^ - A^iA'^) = J"^ . (22) 

This equation describes the possible singularities of the Goldstone field since the first term 
translates as the curl of the gradient of (p. Using the Stokes theorem, the flux of J°^ is 
related with the circulation / V0 ds = kn around the vortex lines, where k is an integer. 
This condition is imposed in the Kalb-Ramond field scheme hj rj = kj , where k is the 
vortex topological number. 

Equation (^) provides the general relation between current and vorticity 

^E^^^-D, (A;,i {r - A'^)) = J"^ , (23) 

and allows us to infer the vortex distribution from the current. Using the Stokes theorem 
the circulation | / dx'^X^^ {J" — A'^) in the perimeter is equal to the flux / dx^^dx'''E^^ai3J°'^ 
over the enclosed surface, and we can write 
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A;,^ {r - A'^) dx^ = 7cN, , (24) 

where the second term is the flux of J""^, and represents the total number of vortices enclosed 
in the integration path. 
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The presence of a vorticity gives an additional term to the force acting on the fluid. The 
force is given by 

r = D,T^^ = -J.^H^^'' - rO^ {\-'^X-) + ^rj-D^ (A;,i) , (25) 

where the first term on the right hand side comes from the new term in the equation of 
motion for H, and the Bianchi identities have been used. The last terms are only relevant in 
regions where the background has a large gradient, with a scale of the order of the coherence 
length, as happens for the dynamics of the vortices in the crust of a neutron star where the 
vortex pins. 

We will now calculate the force acting on a single vortex. In this case the solution of 
(^) could be written as a background solution, whose source are the other vortices, B^'^, 
plus a particular solution with the vortex as a source, -B^'', i.e. H^'^'^ = Hj^^'^ + H!^"'^ . The 
term H^"^ can be written using relation (1T6|), from the average current 

^^au ^ _E^'-'^PJp . (26) 

We are not taking into account the dynamics of the vortex, which could be responsible for 
tension effects but is not given by the effective theory. These effects are related to interaction 
terms between J^^cr and H^^'^, which are not relevant if the curvature radius of the vortex is 
much greater than the coherence length. Considering also that the background potentials 
have a slow variation along the vortex size, the force on the vortex per unit length becomes 

r = -r ^ -j^aE'^'"''.!, . (27) 

In general for a vortex in equilibrium the Magnus force compensates the total external forces 
acting on it. Considering only the interactions here discussed, the equilibrium is specified 
by the null Magnus force condition, which can be put in terms of the current as 



J" 



Ds A;: (J'^ - A-) - D^ (X-^ [r -X'^)) = . (28) 



The vortex can be locally characterized by its tangent vector my and its quadrivelocity 
v'^ . Thus we have 
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da^"" = {v^m" - m^t;^) dadr , (29) 

where a is defined such that m^m^ = 1 and r is the proper time on the vortex, so v^v^, = —1. 
Using J^ = UqU^, where u^ is the fluid quadrivelocity, and integrating over a transversal 
cut, the force per unit length flnally results 

f^ = -^noV^g^%,.rU''v''m^ . (30) 

This is the expression for the generalized Magnus force. Introducing Su^^ = f ^ — -u^ and 
considering that to flrst order in the velocities 6u^ = 0, the 3-d force acting on a vortex with 
m^ = (0, m) is 



TT 



f = 2 ^0 v^V -^°° [G- {6nx m)] , (31) 

where G = {g''^), and f = (/*), and the power dissipated by the interaction of the vortex 
with the fluid is 



f = ^noV^^f^ (V^v.u X m + 2g.(5ux mj . (32) 

In flat space-time this force reduces to the ordinary Magnus force f = |no (5uxm). 

The Magnus force implies that the equilibrium is reached when the vortex moves at 
the same velocity as the medium, or when the vortex is aligned with the relative velocity 
between the vortex and the background. 

IV. EXAMPLES 

In this section we present two examples of superfluids in a gravitational fleld. The flrst 
is a superfluid in a laboratory on Earth. It is simple, but its relevance is only academic 
because the interesting effects are extremely weak. On the other hand, the second example 
is of great phenomenological interest, because it corresponds to a superfluid in a neutron 
star, which is one of its main components. 

In the construction of the effective Lagrangian theory we include all the terms consistent 
with the symmetries of the problem and relevant for the energy scale that matters. In 
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general there is more than one Goldstone boson that arise from the symmetry breaking. In 
the following examples we consider only one boson. This is true for He^ superfluid but not 
for He^ or the neutron star case, where part of the superfluid can be in the spin two ^P2 
state. In spite of this the baryonic phase can be factorized from the spin degrees of freedom. 
Specifically, for the neutron star the effect of the interaction of the baryonic Goldstone boson 
with these spin degrees is negligible ^. This is easy to see by considering that the ratio 
between the spin angular moment and the orbital angular moment for a single neutron is 
^^2^ -C 1, indicating that the spin can not significantly alter the global dynamics of the 
baryonic current. 

A. Superfluid in a laboratory on Earth 

In this case we have a weak gravitational field without a noticeable contribution from 
the fluid itself, with a time independent metric that can be written: 

^°° = 1 - 2f/ , (33) 

g'i = -{! + 2U) 5'^ , (34) 

g'' = h' , (35) 

where U = g z is the Newtonian potential and h gets contributions from the dragging 
and inertial effects due to the Earth rotation. The existence of an energy extremum for a 
null gradient implies that A^ is proportional to the time derivative Killing vector, with 
Ao = n{z) and Aj = 0. Here the minimum energy condition is equivalent to X^^ > and to 
the matrix A*-' negative definite. In general both the gravitational dragging and the fluid 
velocities are small, so that the second order tensor A^'^ can be considered independent of 
h. Taking into account the rotational symmetry on the vertical z axis, the traslational 
symmetry in the horizontal (x, y) plane, and the invariance under temporal inversion, the 
non null components are 

A° = a{z) ,Xl = Xl = Xl = b{z) . (36) 
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Here a and b are scalars that depend only on the fluid conditions. We take \^ = X^ = 
XI because the gravitational fleld scale is much longer than the microphysics scale that 
determines the superfluid state, so that any possible microscopic anisotropy induced by the 
gravitational field should be negligible. The current can be written as 

jo = a<i) + n , ji = bdicf) . (37) 

Let us consider the stationary case. The charge density in the fluid rest frame no is equal 
to n at flrst order in the velocity. The tridimensional velocity is given by the contravariant 
components 

v=—{f) = h- — {l + 2U)V(p. (38) 

rio rio 

This implies V x ( "'°^ ~ — '-(v — h)) = 0, a generalization of the usual superfluid relation 

V X {7= 0. This equation is analogous to the superconductor equation V x (v — ^A) = 0, 
but there is a subtle detail in this analogy. In our case the origin of the velocity drift is 
the metric and not the connection, as in the superconductor case. Thus for a superfluid the 
velocity fleld in the fundamental state will not be null, but equal to h. The term with the 
gradient of the fleld in Eq. ( |58| ) is non null when there are vortices contributing to the fluid 
rotation. As can be deduced from this equation the contribution of h generally diminishes 
the number of vortices with respect to a flat space situation. 

When h has an inertial origin, due to choosing a rotating reference frame, we have 

V X h = —w, the angular velocity vector. This tells us that the superfluid is in fact non 
rotating. Inertial effects on superfluids have been used in the construction of a precision 



gyrometer |11]. In the case of a superconductor, inertial effects lead to a compensating 



London fleld ||T2[, which is also induced by the gravitational dragging ||13|| . 

The effective mass of the boson condensate m* can be identified from Eq. (^) at null 
gravitational field, as the coefficient of V0. It results m* = Uo/b. The Lagrangian density is 

L = {1-2U) ^002 - b (V(j)f + n^] + (2a0 + n) (h ■ V0) - AbU (v^)^ . (39) 
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The limit of zero gravitational field allows us to clarify the physical meaning of n, a and 6, by 
comparing with the low energy BCS Lagrangian |]T^. Thus n = kp/Sir"^, with kp the Fermi 
wave number, is the charge density. We can also relate a and b with the Fermi wave number 
and the fermion mass m, a = ■mkF/'^n'^ and b = kp/Gn'^'m, and thus the sound velocity 
is given by J^ = — j= = vp/^/S, where vp is the Fermi velocity. The relation between 
the fermion mass and the effective pair mass is m* = 2m. Finally, the factor (1 — 2U) 
corresponds to a redshift effect. 

If there are vortices, the expression (^) gives the Magnus force acting upon them. In 
this case we have the usual flat space-time force only affected by a redshift factor. 

B. Superfluid neutron star 

Here we develop with some detail the example of the superfluid component in a neutron 
star. In the case of a neutron star at rest we can assume that the metric satisfies a spherical 
symmetry 

(dsof = e^^dt^ - e^^dr^ - r\de^ + sin ^6 d^^) . (40) 

The potentials $(r) and A(r) that characterize this metric are given by the Oppenheimer- 
Volkoff equations for cold stars [|l^ . 

In the case of neutron stars the friction is high enough to drive the star to equilibrium 
in relatively short times. This allows us to assume a stationary neutron star, which rotates 
at a small angular velocity Vt. This magnitude represents the angular velocity measured by 
an observer at rest with the fluid, and is related with the fluid velocity by fi = m'^'/m*, such 
that 

j'^ = n„n°(l,0,0,(]). (41) 

We can make an expansion of the rotating star metric in the perturbation with respect 
to the rest star metric and the angular velocity Vt. The source for this perturbation is the 
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superfluid energy inomentum tensor. It can be considered a small quantity for ordinary 
pulsars because the gravitational acceleration, approximately ^-, is much greater than the 
centrifugal one, Q^R. As an example, for the rapidly rotating pulsar Crab the quotient of 
these accelerations is approximately 50, implying small deformations. Hence we can keep 
the expansion up to the first order in the fluid angular velocity |]16 



ds'^ = {dsof - 2r2 sin ^9 udipdt. (42) 

The gradients (p^y can also be considered first order quantities because they can be chosen 
to be zero for the star at rest. 

In the construction of the effective Lagrangian theory we will consider all the terms that 
satisfy the symmetries of the problem and are relevant for the energy scale that matters. 
We can identify several symmetries: the baryonic global symmetry, the covariance under 
general coordinate transformations, the spherical symmetry, the temporal translation and 
the temporal reflection symmetry of the star at rest. 

The time translation and rotational invariance makes A^ = /(r) x^, where X/i = 
(e^*, 0,0,0) is the time Killing vector. The covariant components of this Killing vector 
are non dynamical, contrary to the contravariant components, which depend on the metric 
perturbations. 

Time translation and spherical symmetry impose important restrictions to A^*^. We must 
have 

\l = a{r) , a; = b{r) , A^ = c{r) , A^ = c{r) , (43) 

where a, b and c are scalar quantities, because they are the eigenvalues of a mixed tensor. 
The stability of the ground state requires that a, b, and c must be positive. The metric 
perturbations do not affect the quadratic term in o- because this is already of second order. 
The interpretation of b and c comes from the equations 

b(r) 1 

Ur = -^^dr(p = —^fTdr(p , (44) 

no[r) m*[r) 
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U0 = ^de<P = -4r^de<P , (45) 

no(r) mgir) 

C\Tj 1 

u^ = -^^d^(j) = ^—^d^(p , (46) 

nQ[r) m*[r) 

where m* = no/&, Tn*^ = m* = uq/c are scalar quantities of dimension one that depend 
only on the fluid conditions, and because of this are functions of the radius. By using 
the equivalence principle, these parameters can be interpreted as the effective mass of the 
bosonic quasiparticles in the different directions. This is the double of the fermionic particle 
effective mass, which can be calculated from the particular model as the ratio between the 
Fermi momentum and the Fermi velocity. 

The quantities that determine the fluid conditions are related to each other independently 
of the gravitational field. This is because this field can be considered constant along the 
microscopic nuclear interaction range that determines the plasma dynamics. Because of 
this, and as long as the Fermi surface of the fluid has no preferred direction, we can take 
m*j. = m*Q = m*^ = m*. For a neutron star m* includes plasma effects and has values of the 
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For the star at rest we can identify jo = -^o = "^oC*. Up to flrst order in u, for a rotating 
star uq is the same as for a star at rest, and thus this equation is also valid for a rotating 
star up to flrst order in the angular velocity. 

The relation between the angular velocity Q and the metric coefficient u in absence of 
vortices is given by Eq.(^. The solution of this equation is fi = uj{r,6) , i.e. the angular 
velocity of the fiuid must be the same as the angular velocity that takes an object that falls 
free from infinity to the point (r, 9), and corresponds to the angular velocity of local inertial 
frames with respect to the fixed stars. 

Furthermore, in the case we are discussing here, the star that produces the gravitational 
field is formed by the rotating fiuid. Thus u must satisfy the Einstein equation for the 
R^t Ricci tensor component, with the energy momentum component T^t of the fiuid as the 
source. The solution for the perturbed Einstein equation is a 6'-independent u that satisfies 
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UJ,rr + (- - A' - $' ) UJ,r +' (- + $'" A' " -c'^) UJ 

ly ly /i" /y 

= 8ne^^{{p + 3p)u -2{p + p)Q) . (47) 

For nonsingular gravitational potentials A and $ the only solution for Q = u with regular 
geometry is cj = [^. That means that there is no rotating star solutions without vortices, 
i.e. a rotating superfluid star necessarily contains vortices, in which case {Q — uj) is non 
null. The background acts on these vortices with a force given by Eq. (pT]). If this is the 
only force upon the vortex, from it we can obtain the vortex profile in the stationary state. 
In this case the vortex adapts to the background in such a way that this force becomes 
null. This configuration is given by Eq. (|28|) , and in the case we are discussing here it is 
consistent with j^ = A° = noe~'^ and j'^ — X^ = rioC'^ {Q ~ u). Furthermore, Eq. (^) gives 
the relation between the current j^ and the vorticity J'^'^, and allows us to compute the 
distribution and orientation of the vortices. According to Eq. (p^, in this case we have 



n-uir,e)= ^"^''l?, , \, , (48) 

V ' ^ m*e-*Mr2sin2^ ^ ^ 

where Ny is the number of vortices within a closed circular path that passes by the point 
(r, 6) and is perpendicular to the z axis, m* is the effective mass, and e~* gives the red shift 
of the effective mass due to the gravitational field. 

For a rotating star at the minimum energy configuration Q is constant, and thus Eq. 
(^) together with Eq. (^) allows us to completely determine the vortex profile. The angle 
f3 between the axis of the star and the direction of the vortex lines at a given point (r, 6) 
can be deduced from Eq. (^8]), using the cylindrical symmetry of the system. It reads 



. ^ 1 K(r) sin 2^ 

sin/3 = ^^ ^^ -75 , 49 

2(l + ft:(r)(K(r)-2) sin2(^))i/2 ' v J 

where K{r) = | (fXi^^ ~ ^^^ + ^)- '^^^ ^^^^ term in k corresponds to a purely geo- 
metrical contribution due to the Lense-Thirring effect, whereas the two last terms depend 
on the microscopic structure of the fluid and the star as a whole. 

In particular, if we take r equal to the radius R of the star and 6* = |, we have cu = ^^, 



where J is the angular momentum of the star, and thus the total number of vortices contained 
by the star is 

N = m*e-*(«) (qR' - ^) . (50) 

A noticeable effect of the gravitational background is a decrease of the vortex density with 
respect to a similar situation in flat space-time. In the case of a typical neutron star this 
decrease is of the order of 15%. 

V. FINAL REMARKS 

We have considered the superfluidity phenomena in presence of a gravitational back- 
ground. Our starting point is the spontaneous symmetry breaking of a U(l) symmetry. This 
field-theoretical approach allows us to develop from first principles a fully covariant formal- 
ism. Within this framework we study general aspects of the dynamics of the superfluid in a 
gravitational field. We analyze the force acting on a vortex, which is a generalization of the 
well known Magnus force, and in particular we find the profile of a vortex in equilibrium with 
the condensate. This approach makes contact with the microscopic theory which describes 
the details of the superfluidity phenomena in two points. One is a very basic and general 
one, with a strong theoretical support, which is the spontaneous symmetry breaking mech- 
anism. The other point of contact is the specific description of the background through the 
tensors X'^'^ and A^, which in general must satisfy some symmetry requirements, but their 
details depend on the microscopic physics of the system. In some relatively simple cases it is 
possible to construct a rather closed description, dependent only on a few phenomenological 
parameters, making use of the known symmetries. 

Our results are exemplified with two systems. One of them is a superfluid in a terrestrial 
laboratory, where we have a weak gravitational field. The main result here is the formal 
analogy with a superconductor, with the gravimagnetic field h playing the role of the vector 
potential A, and the Newtonian potential U the role of the Coulomb potential. This analogy 
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is not complete, because the geometrical origin of these fields is different. In the gravitational 
case h is introduced by the metric, whereas the A field is due to the connection. One 
consequence of this difference is that in the electromagnetic case we have a Meissner effect, 
while in the gravitational case we have an anti-Meissner one [^. 

The other example considers a more interesting system, which is the superfluid in a neu- 
tron star. In this case we have a strong gravitational field, which makes a fully covariant 
treatment unavoidable. Here we construct a description of the system based on the symme- 
tries of the star, and this is enough to determine the shape and distribution of the vortices, 
assuming that there are no external forces. This should be the case in the star core where 
there are no pinning forces. Other forces, such as the magnetic ones, cannot appreciably 
alter the vortex distribution because their energy is much smaller than the rotational energy. 
The generalized Magnus force here analyzed could be very relevant for a detailed study of 
the transient during the pulsar glitches, if we have an adequate model for the pinning forces. 
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